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SUMMARY 

The application of the finite element method (superelement technique) In 
conjunction with basic concepts from mechanics of materials theory Is demon- 
strated to predict the thermomechanical behavior of high temperature metal 
matrix composites (HTMMC). The simulated behavior Is used as a basis to estab- 
lish characteristic properties of a unidirectional composite Idealized as an 
equivalent homogeneous material. The ply properties predicted Include: ther- 

mal properties (thermal conductivities and thermal expansion coefficients) and 
mechanical properties (moduli and Poisson's ratio). These properties are com- 
pared with those predicted by a simplified, analytical composite micromechanics 
model. This paper Illustrates the predictive capabilities of the finite ele- 
ment method and the simplified model through the simulation of the thermo- 
mechanical behavior of a PlOO-graphI te/copper unidirectional composite at room 
temperature and near matrix melting temperature. The advantage of the finite 
element analysis approach Is Its ability to more precisely represent the com- 
posite local geometry and hence capture the subtle effects that are dependent 
on this. The closed form micromechanics model does a good job at representing 
the average behavior of the constituents to predict composite behavior. 


INTRODUCTION 

High temperature metal matrix composites (HTMMC) are emerging as candi- 
date materials offering high potential payoffs 1n hot structures applications 
for advanced aerospace propulsion systems. Realization of these high potential 
payoffs depend on the development of: (1) fabrication techniques, (2) experi- 

mental techniques for characterizing the complex nonlinear behavior anticipated 
for these applications, and (3) computational methods for engineering design 
that account for this complex behavior. 

Recent research at NASA Lewis Research Center has demonstrated two 
approaches for predicting HTMMC behavior. These approaches consist of: 

(1) application of the finite element method (superelement technique) In con- 
junction with concepts from mechanics of materials theory and (2) application 
of closed-form equations derived from a simplified, analytical micromechanics 
model. The first approach has previously been demonstrated by Caruso (ref. 1) 
and Caruso and Chamis (ref. 2) to predict equivalent thermal and mechanical 
properties of polymer matrix composites. The second approach employs a varia- 
tion of the micromechanics model developed by Hopkins and Chamis (ref. 3) In a 
form more recently Implemented by Hopkins and Murthy (ref. 4). 

The objectives of this paper are to: (1) describe the application of the 

finite element method In conjunction with the mechanics of materials concepts 
as an approach to HTMMC micromechanical modeling and (2) compare the finite 
element results with those predictions based on the closed-fprm model Imple- 
mented In the Metal Matrix Composite ^alyzer (METCAN computer code) (ref. 4). 


LIST OF SYMBOLS 


A area 

dT change In temperature 

E Young's modulus 

F force 

G shear modulus 

K thermal conductivity 

L length 

qt flux 

u displacement 

a thermal expansion coefficient 

V Poisson's ratio 

Subscripts 

x.y.z denotes the components of the global structural Cartesian coordinate 
system 

1,2,3 denotes the components of the composite material Cartesian coordinate 
system 


MODEL DESCRIPTION 

In the analytical micromechanics model, the representation of the funda- 
mental unidirectional composite unit cell Is a square array. The unit cell 
contains a single circular fiber surrounded by matrix to form a square array. 

In the finite element model, the representation Is also a unit cell. The unit 
cell subsequently becomes the "superelement" that can be "Imaged" to form the 
complete nine-cell model (fig. 1). The Impetus for adapting the super- 
element technique Is derived from: (1) the ability to represent a single 
square array (fiber/matrix unit cell) that could be easily modified, (2) the 
consistency with which the single unit cell can serve as the building block to 
form a multicell model with the ability thereby provided to create a stiffness 
matrix of the multi cell model with much fewer total degrees of freedom compared 
to the full discretization, and (3) the capability to simulate multicell models 
substantially larger than conventional finite element analysis would allow. 

The superelement model (fig. 1) consists of one primary superelement (center 
cell) containing 192 three-dimensional solid (hexahedron and pentahedron) Iso- 
parametric elements and eight secondary or "Image" superelements, comprising a 
total of 1133 grid points. The MSC/NASTRAN (ref. 5) general purpose finite 
element analysis software was used for the simulation. The load and boundary 
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conditions are applied to the model as enforced grid point displacements and/or 
temperatures. 

The single cell Is modeled In such a manner so as to allow the fiber 
volume ratio (FVR) to be varied to any one of the four values: 7, 22, 47, or 
62 percent, rather routinely. This Is achieved by changing the elements along 
the circumference of the fiber to matrix elements thereby decreasing the diam- 
eter of the fiber. As previously mentioned, the particular composite system 
chosen for the study consists of PlOO-graphI te fiber (an orthotropic material) 
and a pure copper matrix (an Isotropic material). Values of constituent pro- 
perties at room (70 °F) and elevated (1500 “F) temperature assumed for the 
study are summarized In table I. 


ANALYSIS PROCEDURE 

The following sections describe the analysis procedures to determine the 
mechanical and thermal properties, which Involves the choice of appropriate 
boundary and loading conditions together with the relevant equations from 
mechanics of materials concepts. Composite moduli are determined by obtaining 
forces from the finite element stress analysis resulting from prescribed dis- 
placements and calculating the property using the associated mechanics of mate- 
rials equations. Likewise, composite Poisson's ratios and thermal expansions 
coefficients are calculated by obtaining the displacements from the finite 
element stress analysis which results from prescribed mechanical and thermal 
loading, respectively, again calculating the property using the associated 
mechanics of materials equations. The composite thermal conductivity is 
determined by obtaining the flux and free surface temperature from the finite 
element heat conduction analysis and calculating the composite property from 
thermodynamics concepts. The composite properties at elevated temperatures 
are predicted assuming steady state equilibrium conditions. Under these condi- 
tions, temperature dependent properties and linear analysis are considered to 
be adequate. 


MECHANICAL PROPERTIES 

The mechanical properties simulated from finite element analysis and the 
predictions of METCAN methods are described below. The properties are for 
room temperature (70 “F) and high temperature (1500 °F) conditions. 


Longitudinal Modulus and Poisson's Ratio (Ej^n and 

The boundary and load conditions for the finite element simulation are 
applied to the composite through enforced displacements (Including prescribed 
zero displacements). The boundary conditions are as follows: 

(1) Constrain the degrees of freedom on the front YZ-plane In the 
X-dIrectlon (fig. 2). 

(2) Constrain the degrees of freedom on the centerline parallel to the 
X-axis on the top and bottom XY-planes In the Y-dIrectlon. 
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(3) Constrain the degrees of freedom on the centerline parallel to the 
X-axis on the right and left XZ-planes In the Z-dIrectlon. 

The composite Is loaded with a prescribed uniform displacement on the back 
YZ-plane In the X-directlon. 

The composite longitudinal modulus Is calculated from the following 
equation: 


^ail 


where Is the applied stress due to an enforced displacement on the YZ 

face and Is the longitudinal strain. 

The longitudinal modulus simulated by the finite element analysis Is plot- 
ted as a function of fiber volume ratio (FVR) for the room and high temperature 
conditions In figure 3. The METCAN predictions are plotted for comparison pur- 
poses In figure 3. The comparison between the two results Is, as expected, 
very good. 

The longitudinal Poisson's ratio Is calculated from the following 
equation: 




where Is the longitudinal strain In the load direction and 

the transverse strain In the unloaded direction. 

Poisson's ratio as predicted by the finite element simulation shows little 
effect to temperature (fig. 4). Poisson's ratio Is affected by neighboring 
fibers and the multicell finite element results show this when compared to the 
METCAN predicted result which does not take this affect Into account. The 
affect of other fibers on Poisson's ratio represents a displacement compatibil- 
ity (due to Poisson effects) between the cells which Is not considered for the 
single cell model or METCAN. The METCAN prediction becomes worse as the FVR 
Increases. This Is because as the FVR Increases the displacement Incompatibil- 
ity becomes more predominant; however, the trend Is In the same direction. 


Transverse Modulus and Poisson's Ratio (£ 5,22 va 23 ) 

The boundary and load conditions for the finite element simulation are 
applied through enforced displacements (Including prescribed zero displace- 
ments). The boundary conditions are as follows: 

(1) Constrain the degrees of freedom In the center YZ-plane In the 
X-dIrectlon (fig. 2). 
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(2) Constrain the degrees of freedom In the bottom XY-piane in the 
Z-d1rection. 

(3) Constrain the degrees of freedom on the centerline parallel to the 
X-axis on the top and bottom XY-plane in the Y-direction. 

The composite is loaded with a prescribed displacement on the top XY-plane in 
the Z-direction. 

The transverse modulus is calculated from the following equation: 


F = ^^33 


where the applied stress due to an enforced displacement on the YZ 


face and ^j^33 


is the longitudinal strain. 


As shown in figure 5, the modulus decreases with Increasing FVR. This is 
expected because the transverse modulus of the matrix is larger than that of 
the fiber. The METCAN predictions (the continuous curves) are practically 
identical with the finite element simulations at both temperatures but deviate 
slightly with decreasing FVR. 


The transverse Poisson's ratio is calculated from the following equation: 


_ !M2 
''S123 = c,33 

where aj ^33 is the transverse strain in the load direction and ^^22 
transverse strain in the unloaded direction. 

Poisson's ratio as predicted by the finite element simulation shows little 
effect to temperature (fig. 6). The Poisson's ratio (vj^23^ greatly affected 

by the neighboring fibers and the multicell finite element results show this 
when compared to the METCAN results which do not take this affect into account. 
Values for Vg ^23 plotted assuming transverse isotropy (curves labeled 

isotropic) for comparison to the two methods described here. That is the trans- 

( ^2,22 ,1 
Vg^23 * 2G " ' 1 

not ue jubiiiieu on trie constituent levei. 


Shear Modulus Gai2 

The boundary and load conditions for the finite element simulation are 
applied through enforced displacements (including prescribed zero displace- 
ments). The boundary conditions are as follows: 
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(1) Constrain the degrees of freedom In the left XZ-plane In the 
X-dIrectlon (fig. 2). 

(2) Constrain the degrees of freedom for all nodes In the Y- and 
Z-dIrectlon. 

The composite has a prescribed displacement on the right XZ-plane In the 
X-dIrectlon, 

The modulus Is calculated from the following equation: 

r iBJi 

ill2 Yjii2 

where Is the applied stress due to an enforced displacement on the XZ 

face and Yj^^2 transverse strain. 

The finite element predictions (fig. 7) show the high temperature property 
Is smaller than the room temperature property. This Is expected since at high 
temperature P100 and copper have lower shear moduli than at room temperature. 
The METCAN predictions (continuous curve) compare well with the finite element 
prediction at both temperatures through the full range of FVR. These results 
further Indicate that adjacent fibers have negligible Influence on the shear 
modulus. 


Shear Modulus Gft 23 

The boundary and load conditions for the finite element analysis are 
applied through enforced displacements (Including prescribed zero displace- 
ments). The boundary conditions are as follows: 

(1) Constrain the degrees of freedom In the left XZ-plane In the 
Z-dIrectlon (fig. 2). 

(2) Constrain the degrees of freedom for all nodes In the X- and 
Y-dIrectlon. 

The composite has a prescribed displacement on the right XZ-plane In the 
Z-dIrectlon. 

The modulus Is calculated from the following equation: 

r 

'23 - Y,23 

where x ^22 applied stress due to an enforced displacement on the XZ 

face and Yj^23 transverse strain. 

The finite element predictions for Gg ^23 similar to those for Gj ^^2 

(fig. 8). The METCAN predictions (continuous curve) are In good agreement with 
those from the finite element method at both temperatures through all FVR. 
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THERMAL PROPERTIES 


The thermal properties synthesized from the finite element simulation and 
those predicted by the simplified micromechanics model implemented in METCAN 
are described below. The properties are for both room temperature (70 ’F) and 
high temperature (1500 “F) conditions. The thermal properties of interest 
Include expansion coefficients and conductivities. 


Longitudinal and Transverse Thermal Expansion Coefficient 

The boundary and load conditions for the finite element analysis are 
applied through enforced displacements (including prescribed zero displace- 
ments) and temperatures. The boundary conditions are as follows: 

(1) Constrain the degrees of freedom on the center YZ-plane in the 
X-direction (fig. 2). 

(2) Constrain the degrees of freedom on the centerline parallel to the 
X-axis on the top and bottom XY-plane in the Y-direction. 

(3) Constrain the degrees of freedom on the centerline parallel to the 
X-axis on the right and left XZ-plane in the Z-direction. 

The composite is subject to a uniform temperature applied to the entire 
structure. 

The thermal expansion coefficients are calculated from the following 
equations: 


_ !mi 

“ill 1 ■ VT 

_ !M2 
“222 ■ 7T 

where ^ is the longitudinal strain, c ^22 transverse strain and VT 

is the temperature change. 

The longitudinal expansion as shown in figure 9 decreases rapidly with 
increasing FVR. This large decrease is due primarily to the negative expansion 
coefficient of the P100. The METCAN predictions (continuous curve) are in good 
agreement with and follow the same trend as the finite element predictions 
through all FVR. 

The finite element predictions for the transverse expansion (fig. 10) are 
affected by the neighboring fibers. The expansion coefficient increases 
slightly with increasing FVR up to about 0.4 and then the decrease slightly 
with increasing FVR greater than 0.4. This occurs at high temperature and room 
temperature. The METCAN predictions follow the same general trend as the super- 
element predictions at room and high temperature for all FVR. The room temper- 
ature METCAN prediction (continuous curve) is in better agreement than the high 
temperature prediction although both are considered reasonable based on the 
simplified assumptions in METCAN. 
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Longitudinal Thermal Conductivity 

The boundary and load conditions for the finite element analysis are 
applied through enforced temperatures. The boundary conditions are as follows 

(1) Constrain the temperature of the atmosphere (fig. 2) on the cold face 
(front YZ-plane) to 0 °F. 

(2) Insulate the structure on the four sides parallel to the heat flow. 

(3) Constrain the back YZ-plane to conduct heat to the atmosphere. 

The composite is subject to an increase in temperature to the front YZ-plane. 

The longitudinal thermal conductivity is calculated from the following 
equations: 


, ^t ’ ^ 

■^illl " A . 7T 

where, is the flux through the composite, is the X length, A is the 
cross-sectional area, and 7T is the temperature change. 

The finite element and METCAN predictions of the longitudinal conductivity 
are shown in figure 11. The room and high temperature predictions Increase 
with increasing FVR. The room temperature METCAN predictions are slightly 
lower than the finite element prediction while the high temperature METCAN pre- 
dictions are slightly higher than the superelement predictions. 


Transverse Thermal Conductivity 

The boundary and load conditions for the finite element analysis are 
applied through enforced temperatures. The boundary conditions are as follows: 

<1) Constrain the temperature of the atmosphere (fig. 2) on the cold face 
(right XZ-plane) to 0 “F. 

(2) Insulate the structure on the four sides parallel to the heat flow. 

(3) Constrain the left XZ-plane to conduct heat to the atmosphere. 

The composite is subject to an Increase in temperature to the right YZ-plane. 

The transverse thermal conductivity is calculated from the following 
equations: 


S22 



where, q^ is the flux through the composite, is the Y width, A is the 
cross-sectional area, and 7T is the temperature change. 
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The finite element and METCAN predictions are shown In figure 12 for the 
transverse thermal conductivity. The predictions decrease with Increasing FVR 
for room and high temperature. The METCAN results are lower than the finite 
element predictions for both room and high temperature. Since the finite ele- 
ment predictions are higher at low FVR and approach the METCAN results as the 
FVR Is Increased, It may be concluded that denser whiskers may be required for 
smaller diameter fibers. 


SUMMARY OF RESULTS 

The Important results of adapting mechanics of materials concepts with 
finite element analysis, or alternatively the simplified analytical models such 
as Is Implemented In METCAN, to HTMMC micromechanics are as follows: 

(1) The mechanics of materials concepts and the finite element (super- 
element) method provides an independent method of predicting the thermomechani- 
cal properties of MMC at room and high temperature. 

(2) The mechanical properties Include moduli and Poisson's ratio. Compar- 
isons with the simplified micromechanics equations (METCAN) are in good agree- 
ment. Poisson's ratio (vg,,) does not satisfy the assumption of transverse 
Isotropy. 

(3) The thermal properties include thermal expansion coefficients and 
thermal heat conductivities. The longitudinal thermal expansion coefficient 
predicted by METCAN Is In good agreement with the finite element results. The 
transverse thermal expansion coefficient predicted by METCAN deviates slightly 
from the finite element results. The thermal conductivities predicted from 
METCAN are In good agreement with the finite element simulations. 


CONCLUSION 

The Important conclusions of adapting the superelement method to thermo- 
mechanical behavior of HTMMC are as follows: 

1. The finite element technique used in conjunction with mechanics of 
materials concepts provides the ability to analyze ordinary composite struc- 
tures well and offers an alternative for predicting composite properties. 
Because of Its utility this method can be extended to analyze more complex 
composite Interactions. These may include fiber fracture, fiber/matrix dis- 
bonding, hybrid composites, etc. 

2. The superelement method can provide detailed information about local 
composite behavior that cannot be obtained through experiments. 
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TABLE I. - FIBER AND MATRIX PROPERTIES 


Property symbol 

Uni ts 

Temperature, 

°F 



70 

1500 



P100 

Copper 

PlOO 

Copper 

Longi tudi nal 

Mpsi 

105.0 

17.0 

93.8 

8.87 

modul us 





Transverse modulus 

E 22 Hpsi 

.90 

17.0 

.804 

8.87 

Longitudinal shear 

G ]2 ^psi 

1.10 

6.54 

,982 

3.41 

mosulus 

Transverse shear 

.70 

6.45 

.625 

3.41 

G 23 Mpsi 

modulus 
Longi tudi nal 

.200 

.300 

.179 

.300 

V 12 

Poisson's ratio 
T ransverse 

.250 

,300 

.223 

.300 

'^23 

Poisson's ratio 





Longitudinal heat 
conduction 

Ki, Ptu in, 
hr in . 2 

25. 

19.3 

24.4 

16.8 

Transverse heat 

1 ^-- Btu in. 





conduction 

hr °F in. 2 

1.74 

19.3 

1.70 

16.8 

Longi tudi nal 
thermal expansion 

“11 xio"® . 

' ' 1 n. ®F 

-.90 

9.8 

- 1.01 

19.6 

Transverse thermal 
expansion 

“?? -int 

in. ®F 

5.6 

9.8 

6.27 

19.6 
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FIGURE 1. - 9-CELL FINITE ELEMENT MODEL. 
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P100-6RAPHITE FIBER/COPPER MATRIX 


p CENTER LINE 
i OF TOP 



ftTHOD TEMPERATURE, 

°F 

A 3-D FINITE ELEMENT 1500 

METCAN 1500 

□ 3-D FINITE ELEMENT 70 

METCAN 70 



0 


.2 .4 .6 

FIBER VOLUME RATIO 


.8 


FIGURE 2. - GENERAL BOUNDARY CONDITIONS FOR THE 3-D SUPER- 
ELEMENT MODEL. 


FIGURE 3. “ LONGITUDINAL MODULUS OF UNIDIRECTIONAL 
COMPOSITES. 
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P100-GRAPHITE FIBER/COPPER MATRIX 


P100-GRAPHITE FIBER/COPPER MATRIX 



METHOD 

TEMPERATURE 



°F 

A 

5-D FINITE ELEMENT 

1500 


- METCAN 

1500 

□ 

3-D FINITE ELEMENT 

70 

— — 

METCAN 

70 




METHOD 

TEMPERATURE 



°F 

A 

3-D FINITE ELEMENT 

1500 


• METCAN 

1500 

□ 

3-D FINITE ELEMENT 

70 

— 

‘ METCAN 

70 


□ 3 



I 1 ^ ^ J 

0 .2 .4 .6 .8 

FIBER VOLUME RATIO 


0 .2 .6 .8 

FIBER VOLUME RATIO 


FIGURE 4. - POISSONS RATIO OF UNIDIRECTIONAL COMPOSITES. 


FIGURE 5. - TRANSVERSE MODULUS OF UNIDIRECTIONAL COM- 
POSITES. 
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P100-GRAPHITE FIBER/COPPER MATRIX 



METHOD 

TEMPERATURE 



°F 

A 

3-D FINITE ELEMENT 

1500 


METCAN 

1500 


3-D FINITE ELEMENT, 

1500 


ISOTROPIC 


O 

3-D FINITE ELEMENT 

70 

— 

f^TCAN 

70 



3-D FINITE ELEMENT, 

70 


ISOTROPIC 



PI 00-GRAPHITE FIBER/COPPER MATRIX 


METHOD 

TEMPERATURE 



3-D FINITE ELEMENT 

1500 

ftTCAN 

1500 

3-D FINITE ELEMENT 

70 

METCAN 

70 


3 


CO 

k 


CO 


Q 





0 


.2 .il .6 

FIBER VOLUft RATIO 


.8 


FIGURE 6. - POISSONS RATIO OF UNIDIRECTIONAL COMPOSITES. FIGURE 7. - SHEAR MODULUS OF UNIDIRECTIONAL COMPOSITES. 
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MODULUS (6,,,) MPSl 



0 .2 .6 -8 
FIBER VOLUME RATIO 


FIGURE 8. - SHEAR MODULUS OF UNIDIRECTIONAL COMPOSITES. 



0 .2 -6 .8 
FIBER VOLUME RATIO 


FIGURE 3. - LONGITUDINAL THERMAL EXPANSION OF UNIDIRECT- 
IONAL COMPOSITES. 
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FIGURE TO. - TRANSVERSE THERMAL EXPANSION OF UNIDIRECT- 
IONAL COMPOSITES. 


P100-GRAPHITE FIBER/COPPER MATRIX 


FffflOD TEMPERATURE, 

°F 

A 3-D FINITE ELEMENT 1500 


METCAN 1500 

□ 3-D FINITE ELEftNT 70 



0 .2 A .6 ,8 

FIBER VOLUME RATIO 

FIGURE 11. - LONGITUDINAL THERMAL CONDUCTIVITY OF UNI- 
DIRECTIONAL COMPOSITES. 
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CONDUCTIVITY (K, 


P100 GRAPHITE FIBER/COPPER MATRIX 


METHOD 

TEMPERATURE 


°F 

3-D FINITE ELEMENT 

1500 

METCAN 

1500 

3-D FINITE ELEMENT 

70 

METCAN 

70 


3 



0 .2 .If .6 .8 

FIBER VOLUME RATIO 

FIGURE 12. - TRANSVERSE THERMAL CONDUCTIVITY OF UNI- 
DIRECTIONAL COMPOSITES. 
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